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rABLE OF SYMBOLS 





2 




A 


Feet 


Area of foil 


F 


2 




A 

S 


Feet 


Area of strut 


CG 




Center of gravity 


C 


Dimensionless 


Drag coefficient of a foil 


DF 




C 


Dimensionless 


Drag coefficient of a strut 


DS 




C 


Dimensionless 


Lift coefficient 


L 






C 

S 


Dimensionless 


Side force coefficient 


F 


Pounds 


Craft 1 s centrifugal force 


CENT 




F 


Pounds 


Drag force on a foil (water axes) 


DiF 




F 


Pounds 


Drag force on a strut (water axes) 


DiS 




F 


Pounds 


Lift force on a foil (water axes) 


LiF 




F 


Pounds 


Perpendicular force to craft's 


NET 








deck 


F 


Pounds 


Side force on a strut (water axes) 


SiS 




F 


Pounds 


Force in direction of body X-axis 


X 




F 


Pounds 


Force on a foil in body X-axis 


XiF 




F 

Y 


Pounds 


Force in direction of body Y-axis 


F 


Pounds 


Force on a foil in body Y-axis 


YiP 




F 


Poun ds 


Force on a strut in body Y-axis 


YiS 
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F 

Z 


Pounds 


F 

ZiF 


Pounds 


F 

ZiS 


Pounds 


9 x 


Feet/Second 


9 * 


Feet/Second 


9 Z 


Feet/Second 


Helmr 


Radians 


h 

S 


Feet 


I 

XX 


2 

Slug-Foot 


I 

YY 


2 

Slug-Foot 


I 

ZZ 


2 

Slug-Foot 


L 


Foot- Pound 


L 

XiF 


Feet 


L 

YiF 


Feet 



Force in direction of body Z-axis 
Force on a foil in body Z-axis 
Force on a strut in body Z-axis 

Component of gravity in body 
X-axis 

Component of gravity in body 
Y-axis 

component of gravity in body 

Z-axis 
Helm angle 

Height of Height Sensor above 
instantaneous water surface 
Moment of inertia about body 
X-axis 

Moment of inertia about body 
Y-axis 

Moment of inertia about body 



Z-axis 

Moment acting to produce roll 

about body X-axis 

Distance from CG. along body 

X-axis to a point of application 
of foil force 

Distance from CG. along body 

Y-axis to a point of application 
of foil force 
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L 

ZiF 


Feet 


Distance from CG. along body 

Z-axis to a point of application 
of foil force 


L 

ZiS 


Feet 


Distance from CG . along body 

Z-axis to a point of application 
of strut force 


M 


Foot-Pound 


Moment acting to produce pitch 
about body Y-axis 


D1 


Slugs 


Mass of the craft 


N 


Foot-Pound 


Moment acting to produce yaw about 
body Z-axis 


P 


Radians/Second 


Roll rate about body X-axis 


Q 


Radians/Second 


Pitch rate about body Y-axis 


R 


Radians/Second 


Yaw rate about body Z-axis 


Rgir 


Feet 


Craft's turning radius 


Rudr 


Radians 


Rudder angle 


S 

EiF 


Feet 


Submergence of a foil in earth 
axes 


S 

EiS 


Feet 


Submergence in earth axes of the 

assumed point of application of 
strut force 


T 

X 


Pounds 


Thrust in direction of body X-axis 


U 


Feet/Second 


Velocity in direction of body 
X-axis 


u 

E 


Feet/Second 


Velocity along earth X-axis 


U 

i 


Feet/Second 


Foil velocity in direction of body 
X-axis 


V 


Feet/Second 


Velocity in direction of body 
Y-axis 


V 

E 


Feet/Second 


Velocity along earth Y-axis 


V _ 

X 


Feet/Second 


Foil velocity in direction of 
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w 



H 

E 

V. 

i 



X 

E 



Y 

E 



Z 

E 



a 

i 

/? 

i 

p 

* 

e 



Feet/Second 


body Y-axis 

Velocity in direction of body 


Feet/Second 


Z-axis 

Velocity along earth Z-axis 


Feet/Second 


Foil velocity in direction of body 


Feet 


Z-axis 

Distance along earth X-axis from 


Feet 


reference origin 

Distance along earth Y-axis from 


Feet 


reference origin 

Distance along earth Z-axis from 


Radians 


reference origin 

Angle of attack of a foil 


Radians 


Angle of side slip of a strut 


3 

Slugs/Foot 


Water density 


Radians 


Roll angle 


Radians 


Yaw angle 


Radians 


Pitch angle 



S ub scr ip ts 

C Refers to center foil or strut 

P Refers to port foil or strut 

S Refers to starboard foil or strut 

M Refers to raid foil section 

F Refers to foil 

S Refers to strut 

i When i appears in the subscript to a symbol, it 

indicates that the symbol is to be repeated with i 
succesfully replaced by C,P,S and M, as indicated. 
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I* INTRODUCTION 



This thesis deals with the yaw motion and turning 
characteristics, in calm water, of the hydrofoil craft 
PC(H)-1 High Point. This will be the model used in the 
computer simulation program for this study. 

In previous theses, the system was made conditionally 
stable in two modes, the pitch-heave mode and the roll mode. 
In this thesis, keeping the system stable and working in six 
degrees of freedom, its behavior will be analyzed in three 
modes, pitch mode, roll mode and yaw mode, when the 
simulated model is turning by the effects of a helm step 
command, acting on the astern flaps. This analysis will 
also include the study of rudder effects, aft starboard and 
aft port flaps effects and the heeling behavior of the 
system. In a further investigation the forward flap will be 
split in two parts and pitching and rolling effects will be 
analyzed . 

Since slight changes in the manual control of the 
system occur in such short time intervals and provoke 
tremendous effects, the system necessitates an automatic 
control that will provide the necessary outputs to operate 
the model according to the given inputs. This automatic 
control, already existent, will be improved and adapted to 
the new motion of the craft. 

Finally the favourable roll angle for a certain speed 
and a certain helm step command will be calculated and 
compared with the actual roll angle of the craft. 



12 



II* hydrofoil equations s tudy 



A. FOIL TYPES AND FOIL ARRANGEMENT 

There are two principal foil types used in the 
construction of hydrofoil crafts, the surf ace-piercing foil 
and the fully-submerged foil, as shown in Figure 1. 







Surface-Piercing Foil 







FIGURE 1 

TYPE OF HYDROFOIL CRAFT 
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The surface-piercing foil provides a self-stabilization 
for the craft in pitch, height and roll. On the other hand 
the fully-submerged foil cannot adjust itself to a changing 
lift. This can only be achieved by a change in incidence 
angle of the foil flaps. The fully-submerged foil requires a 
more sophisticated sensing arrangement to control that 
incidence. But the fully-submerged foil has a great 
advantage, when the craft is operating in heavy sea. The 
surface-piercing foil responds to relatively small changes 
in the surface of the sea and gives a very rough ride. On the 
other hand, the fully-submerged foil is much less affected 
by the effects of the surface of the sea. This paper will 
only deal with the fully-submerged foil craft. 

The model for this simulation, the High Point 
PC(H)-1,uses the foils in canard arrangement, which has a 
secondary foil as forward foil and a main foil as aft foil. 
Both foils are supported by vertical struts. Figure 2, 
important for the understanding of the digital computer 
simulation program, shows the configuration of the forward 
and aft foil. The drawing and the dimensions of the craft and 
the foils are shown in Figures 3A and 3B. 
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FIGURE 2 

FOILS CONFIGURATION 
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36 



PRINCIPAL DIMENSIONS OF PC(H) 



r’oKr'AKS f&L 





Vl 







$ 
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FIGURE 3B 

DRAWING OF FOILS AND FLAP CONFIGURATION 



B. COORDINATE SYSTEMS 



Three coordinate systems will be used in this study of 
the hydrofoil craft dynamics. They are : the Body Axes, the 
Earth Axes and the Water Axes. Each of these systems is a 
right hand orthogonal system. 

The body axes coordinate system (X,Y,Z), shown in 
Figure 4, has its origin at the center of gravity of the 
craft and is fixed relative to the craft with the X axis 
running longitudinally through the craft. This system is 
initially coincident with the origin of the earth axes 
system, but it moves with the craft as time progresses. X is 
positive forward, Y is positive to starboard and Z is 
positive downward. ■ 




FIGURE 4 



ORIENTATION OF THE BODY AXES WITH RESPECT TO THE 



CRAFT AND DIRECTIONS OF THE POSITIVE VELOCITIES 
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The earth axes coordinate system (X , Y , Z ) is fixed 

E E E 

relative to the earth's surface. The origin may be at any 

convenient location, provided the body system origin is 

initially at the same point. The X axis lies in the 

E 

horizontal plane ; it is initially coincident with, and 
positive in the same direction as, the body X axis. The Y 

E 

lies in the horizontal plane and is positive to starboard 

when the observer is facing the positive X direction. The 

E 

Z axis is normal to the horizontal plane and is positive 
E 

downward. The craft's motion relative to the calm water 
surface, (constraint of this paper) , is described in terms 
of this coordinate system. 

In general, the body axes are displaced from the fixed 
earth axes by the three Euler angles, 0 (pitch angle) , <t> 

(roll angle) and 'F (yaw angle) . Figure 6 defines the 
positive directions of the Euler angles. The transformations 
between earth and body axes are outlined in Ref. 3, by 
transformation matrices, however, only the results will be 
considered here. 

The water axes coordinate system is aligned with the 
relative water velocity vector and resolves into lift, drag 
and side force directions. All hydrodynamic model test data 
is taken and presented as plots in terms of this coordinate 
system. 

The hydrodynamic force and moment data are obtained in 
a watep axis system which is always oriented with respect to 
the relative water velocity. These data must be transformed 
into the body axes in order that their effect will be 
properly included in the equations of craft motions. 

In general, the water axis orientation is always 
changing with respect to the body axes. This leads to the 
need for a transformation matrix, also outlined in Ref. 3, to 
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resolve the lift, drag and side force values in 
axes. The orthogonal water axis coordinate system 
in Figure 5. The angles of rotation between water 
axes are a and /3 and will be defined later. 



the body 
is defined 
and body 



L 




FIGURE 5 

ORIENTATION OF WATER AXES COORDINATE SYSTEM 
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FIGURE 6 

DEFINITION OF POSITIVE DIRECTIONS OF EULER ANGLE 



PI 



C. EQOATIONS OF MOTION 



The development of the six equations of motion for a 
hydrofoil follows the same procedure as those for a 
displacement type hull. This development is clearly outlined 
in Ref. 3 and therefore will not be done here. 

- The equations of motion are summarized in this form : 

Translat ion 



m 

0 

V 

w 



(1/m) F - QW + RV 
X 

(1/m) F^ + PW - RU 

(1/m) F + QU - PV 
Z 



Rot atio n 



m 

P 

m 

Q 

m 

R 



1/1 |L-QR(I -I 

XX l ZZ YY 



) + (R + QP) I 1 
XZJ 



2 2 

1/1 [m-RP(I -I )-(P -R )I 1 

— - L xx zz xzJ 



YY 






( 2 - 1 ) 

( 2 - 2 ) 

(2-3) 

(2-4) 

(2-5) 

(2-6) 



Investigations with the craft simulation show that the 

terms RV,PV, and PW are very small and may be neglected. 

Similarly, the term I is less than 10% of I and less 

XZ XX 

than 4% of I . Consequently, the terms containing I may 
ZZ XZ 

be neglected for this craft. 

Equations (2-1) through (2-6) are linear and rotational 

accelerations in body axes. 

Given the previous equations, the velocity terms in 
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body axes are obtained by integrating each equation one 
time. 



U 

V 

W 

p 

Q 

R 



-/* 



dt 



= /'" 

= /r 



V dt 



dt 



P dt 



Q dt 



R dt 



(2-7) 

( 2 - 8 ) 

(2-9) 

( 2 - 10 ) 

( 2 - 11 ) 

( 2 - 12 ) 



In order to calculate the craft's position in earth 
axes coordinates, these velocity outputs must be transformed 
from body axes coordinates into earth axes coordinates. The 
transformation of the linear equations become : 



U = Ucos 0cos V (cosi/'sinflsin <j> -sin \p cos<i>) 

E 

+ W (cos'/'sin0cos<£+sin'/'sin <t> ) (2-13) 



V = Ucos0sin'/-+ V (cos'/'cos^+sin ^ sin0sin</>) 

E 

+ W (sin'/'sintfcos^-cos sin<£) (2-14) 



W = -Usin0 + Vcos0sin</>+Wcos0cos <j> (2-15) 

E 



Only the craft's velocity in the downward dir ection,Z^, 



is required to calculate its height above the water surface. 
That is, the vertical position in earth axes follows: 

Z = f W dt (2-16) 

E J E 



23 



The craft's position in earth axes coordinates 
obtained from : 

X = f U dt 
E J E 



Y = / V dt 
E J E 



is 



(2-17) 



(2-18) 



The Euler angles, <t> , 9 and are also of interest 

0 

and can be found by integration of Euler rates <f> , 0 and 'P . 
The rates in turn must be derived from the body angular 
rates. The Euler rates are not easily obtained because they 
occur about axes which are not orthogonal. This fact can be 
appreciated by referring to Figure 6 and recalling how the 
angles were defined. '1 is a rotation about the Z axis, 0 
is a rotation about the Y axis, and '$> is a rotation about 
the X axis. They result in the following set of equations : 

<t> = P + iAsin 8 (2-19) 

6 = Qcos4>-Rsin ^ (2-20) 

'P = (Qsin<f> + Rcos</>) cos0+ (<h-P) sinfl (2-21) 



The Euler angles are obtained by integration of these 
equations. 

Studies of this craft indicated that during abnormal 
situations, roll angles of 15 degrees and larger may occur, 
but under normal operating conditions, roll angles will not 
exceed 10 degrees and pitch angles are not expected to reach 
3 degrees. Therefore, the following assumptions are 
justified : 

sin 8 =0 (2-22) 

cos 0 = 1 (2-23 ) 
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D. CALCULATION OF VELOCITY COMPONENTS 



The various foil velocity components must be determined 
at this point in order to calculate the angles of attack and 
side slip. The total velocity components of each foil in 
body axes can be expressed in terras of the craft linear and 
angular velocities by considering the craft geometry as 
represented in Figure 7 and Table I. 

Center Foil 



U = U + L . Q 
C ZCF 



V = V - L .P + L . R 
C ZCF XCF 



W = W - L .Q 
C XCF 



(2-24 ) 
(2-25) 
(2-26) 



Port Foil 



u = u 




L . Q 


- 


L . R 


(2-27) 


p 




ZPF 




YPF 


V = V 


_ 


L .P 




L . R 


(2-28) 


p 




ZPF 




XPF 


w = w 


+ 


L . P 


— 


L . Q 


(2-23) 


p 




YPF 




XPF 


Starboard Foil 












U = U 




L . Q 


_ 


L . R 


(2-30) 


S 




ZSF 




YSF 




V = V 


_ 


L . P 




L . R 


(2-31) 


S 




ZSF 




XSF 




w = w 


+ 


L .P 




L . Q 


(2-32) 


s 




YSF 




XSF 


Mid-Foil Segment 












0 = U 




L . Q 






(2-33) 


M 




ZMF 








V = V 


_ 


L .P 


+ 


L . R 


(2-34) 


M 




ZMF 




XMF 
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W = 


W - 


• L . Q 






M 




XHF 


0=Point of 


application 


of 


foil forces 


0=Point of 


applica 


tion 


of 


strut forces 


L 


= L 








XCF 


XCS 








L 


= L 


L 


= 


II 


XPF 


XPS 


XSF 




XSS XHF 


L 


= L 


L 


= 


L 


YSF 


YSS 


YPF 




YPS 


L 


= L 


L 






ZPF 


ZSF 


ZMF 






L 


and L 


are 


variables and in 


ZPS 


ZSS 









(2-35) 



equal. 







FIGURE 7 

DEFINITION OF SYMBOLOGY 
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Point of Symbol for Force Along Coordinate From CG to Point 

Application of or Parallel to of Application Relative to 
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Thrust 



E. ANGLE OF ATTACK AND ANGLE OF SIDE SLIP 



The angle of attack of a foil is shown in Figure 8A and 
is defined as the angle whose tangent is the total relative 
velocity between the foil and water in the body Z direction, 
divided by the total relative velocity between the foil and 
water in the body X direction. 

o = arc tan (W /U ) (2-36) 

r r 



If the foil is not aligned with the body X axis, then 
the foil has a fixed incidence angle and the total angle of 
attack becomes ; 

a = a + a (2-37) 

total fixed 



The side slip angle is shown in Figure 8B and can be 
represented by : 



0 = arc tan 

i 



, / 2 2 2 ~' 
v An u +v +w ) 

ri » ri ri ri 



(2-38) 



During actual craft operations the angles will be less 
than 10 degrees, so that the inverse trignometric functions 
can be replaced by quotients. Except at low craft speeds, 
the forward velocity of the craft, 0, is much greater than 
the velocities, V and W. Applying these simplifications, the 
equations become : 

a = (W /U) + cl 

i ri i-fixed 

P. = V ./ U 
l r i 
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F. EXPANSION OF FORCES AND MOMENT TERMS 

With the aid of Figure 7 and Table I,the forces and 
moments for the craft are easily written. It is important to 
note that when a numerical value is inserted for one of the 
alphabetical dimensions, the sign of the value must be 
positive or negative according to positive and negative 
directions from the axis system origin at the C.G. One of 
the assumptions made upon commencing simulation was that the 
struts contributed with zero lift and that the foils 
contributed with zero side force. Also note that horizontal 
foils make no contribution to the body Y force and vertical 
struts do not contribute to the body Z force. Therefore, 
all body Y motion originates at the struts and all body Z 
motion originates at the foils. 

The forces and moment equations become : 

F =F +F +F +F +F +F +F +mg +T (2-39) 

X XCF XPF XSF XMF XCS XPS XSS X X 

F =F +F +F + mg (2-40) 

Y YCS Y PS YSS Y 

F =F +F +F + F + mg (2-41) 

Z ZCF ZPF ZSF ZMF Z 

L=(F L ) + (F L ) - (F L ) - (F L ) - (F L ) (2-42) 
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associated with the body axis and requires no further 
expansion. The gravity terms can be simplified using small 
angle approximations for the Euler angles. 

The standard expression for hydrodynamic forces in 
water axis coordinates is : 

F = gAC i=L,D,S (2-45) 

i i 



Where 1,0,3 represent lift, drag and side force respectively. 
The lift, drag and side force coefficients vary as functions 
of angle of attack, angle of side slip, submergence, 
velocity, flap deflections, elevator deflections and rudder 
deflections . 

Ideally, mathematical expressions would have been 
developed to correctly depict the interrelation of all 
variables which affect the force coefficients. However, 
except for the conditions when relationships are linear, 
hydrofoil technology has not advanced to a position which 
would yield such expressions. The difficulties which are 
encountered iji deriving a mathematical expression include : 

1) The ocurrence of cavitation and ventilation, which 
are not completely predictable. 

2) The even more unpredictable cessation of cavitation 
and ventilation. 

3) The nonlinearity of the hydrodynamic coefficients. 

4) Lack of sufficient test data to completely describe 
all of the above. 
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Ill* HYDROFOIL SIMULATION 



A. SIMULATION OBJECTIVES 



The 


objective of this 


simulation 


was 


to study 


the 


dynamic 


characteristics of 


motion along 


and 


about each 


body 


axis. T 


hat is, perturbations should 


be 


introduced 


to 



separately excite : 

1) Motion along the X-axis 

2) Motion about the X-axis 

3) Motion along the Y-axis 

4) Motion about the Y-axis 

5) Motion along the Z-axis 

6) Motion about the Z-axis 

The following assumptions were made prior to beginning 
simulation : 

1) Craft equations of motion are valid only for the 
foilborne mode. 

2) Height of the craft remains constant. 

3) Hydrodynamic coefficients are based on fully wetted 
surfaces, i.e. no cavitation or ventilation. 

4) The craft may be operated in calm water or in a 
seaway. 



32 



B. COMPUTES SIMULATION 



The computer program used was based on the computer 
programs of Sef . 1 and Ref. 2. 

As in the preceding References, before the main program 
could be assembled, Fortran subroutines were used tc find 
the hydrodynamic coefficients of foil lift and drag, and 
strut drag and side slip. The SETUP subroutine stored as 
data points the simulation curves of Ref. 4. The subroutines 
INTERP and INTRP1 are interpolation subroutines to obtain 
the proper coefficient from the above simulation curves. For 
INTERP, values of angle of attack and forward flap angle 
were used to enter the curves and obtain a value of foil 
drag. Angle of side slip and submergence were used to enter 
INTRP1 to obtain strut drag. Listings of the subroutines are 
shown in the computer program section. The main program was 
then assembled using Figure 9 as a reference for data flow 
and is also listed in the computer program section. In one 
of the parts of this main program, the automatic control 
block, were introduced some modifications, needed for the 
main objective of this paper, the study of the craft's 
turning characteristics. This automatic control block is 
still open for future studies. 
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equation relationships 




C. AUTOMATIC CONTROL 

Control signals originate at the pilot's controls, at 
the motion sensors and at the position transducer. These are 
combined and processed in the electronics to produce the 
signal for the servo-valve, The pilot will be able to move 
the helm, a lever, or a knob to introduce steering, altitude 
and attitude commands ; the motion sensors sense errors 
between the commanded and actual craft motions and produce 
electrical signals proportional to these errors ; the 
position transducer provides the feedback signal so that the 
control surface stops moving when it has reached the 
position corresponding to the processed error signals. 

It is now necessary to consider how the control system 
components affect the eguations of motion through the 
control surface deflections. In general, a signal originated 
at the pilot or at the motion sensors will be modified by 
the individual dynamics of the components in the signal path 
before the control surface motion actually occurs. Motion 
sensors to be discussed are vertical gyros, rate gyros, 
accelerometers and height sensors. 

Vert ical Gyros 

A vertical gyro here will produce one signal 
directly proportional to the pitch angle of the craft and 
one signal directly proportional to the roll angle of the 
craft . 

2. Rat e Gyros 

A rate gyro will sense an angular rate about one of 
the three body axes. The dynamics of the gyro are such as to 
produce an output signal which is related to the actual 
angular rate by a second order differential eguation. 

3 • Ac ce l erome ters 

Outputs of accelerometers are also assumed to relate 
the craft motions by second order differential equations. 

To identify the accelerations sensed by the 
instruments, it will be assumed that there is one 
accelerometer at the C. G. , mounted such as to measure 
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lateral acceleration, and one accelerometer above each foil 

mounted in the X-Y plane so as to measure accelerations 

parallel to the 2-axis. 

The lateral acceleration at the C.G. is : 

a = F /m (3-1) 

Y Y 



Taking angular accelerations into account, the 
accelerations in the Z-direction at points above the center, 
port and starboard foil are : 



a = (F /in) -L . Q 
ZC Z XCF 



(3-2) 



a 

ZP 



(F /m) +L .Q-L 
Z XPF YPF 



(3-3) 



a 

ZS 



(F /m) +L .Q+L 
Z XSF YSF 



(3-4) 



4* He ig ht Sen sor 

Ultrasonic and sonar types are either in use or 
expected to be used. For this study it will be assumed that 
the dynamic characteristics of the sensor, whatever the 
type, will be such as to give a direct indication of the 
instantaneous height of the sensor above the water. 
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FIGURE 10 

SPAIF RUDDER CHARACTERISTICS 
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IV* HYDROFOIL TURNING CHARACTERISTICS STUDY 



A. RUDDER EFFECTS 

For the hydrofoil steering control, the struts can be 
viewed as vertical foils so that steering can be achieved by 
rotating a strut, by rotating a hinged trailing edge or by 
rotating an extension of a strut below a foil, (spade or 
ventral rudder), as in the case of the PC(H)-1. 

The rudder on the forward strut exerts an influence 
upon the forward strut side force coefficient. Based on 
Figure 10 the necessary changes were made in the computer 
simulation program. To document the rudder's influence in 
the system and to show the craft's behavior, several figures 
are presented, corresponding to runs using different speeds 
and successive rudder angles of 30,25,20,15,10, and 5 
degrees at each speed. The craft started moving straight 
ahead, the rudder angle being introduced after 3 seconds, 
for all the cases. 

1 . Speed =3 0 Knots 

For a forward speed U = 30 knots. Figure 1 1 A shows tne 
pitch angle of the craft, which had an oscillation at the 
beginning of the run, but it went to a steady state value of 
1.53 degrees, with the craft's bow up, in 4.5 seconds, for 
all the six different rudder angles used. The result was 
natural, because the rudder plays a negligible role in 
changing the value of the foil lift coefficient and 
therefore in the craft's pitch motion. 

For the same speed, in Figure 1 1 B the different roll 
angles correspon ding to different rudder angles were 
represented. When the ship went straight ahead, naturally 
the roll angle was 0 degrees, but as soon as a rudder a ~le 
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of 30 degrees was introduced after 3 seconds, the craft 
started rolling, going to a steady state value of -1.29 
degrees for the roll angle in 19.5 seconds. On the other end 
of the scale for a rudder angle of only 5 degrees, the 
reaction was smaller, the initial oscillation was also 
smaller and the craft went to a steady state value of -0.216 
degrees for the roll angle in 13 seconds. Comparing with 
Figure 4, these negative values show a roll to the port side 
direction, although the ship was turning to starboard. This 
result would be perfectly acceptable for a conventional 
ship, but banking in the opposite direction is desired for 
the hydrofoil because of the high speeds of operation. 

In Figure 11C were represented the several heading 
angles, with positive values, showing that the craft was in 
fact turning to starboard. At 3 seconds when 30 degrees 
rudder angle were introduced, the craft had a reaction, 
reaching a maximum heading angle of 2.02 degrees, but 
returning to a steady heading value of 1.39 degrees in 23 
seconds. For a rudder angle of only 5 degrees the maximum 
heading of the craft was 0.33 degrees , getting a steady 
heading of 0.23 degrees in 17.5 seconds, because the rudder 
change was smaller and then the reaction was smaller. 

• The craft's height above the water Z was 

El 

represented in Figure 1 1 D . Starting with a value of 0 feet, 
the craft went to 3.9 feet above the water in 4.5 seconds. 
Comparing with the dimensions in Figure 3A, it was observed 
that the foils were well inside the water. Again it was 
observed that the rudder plays a negligible role on the lift 
motion of the craft, because the craft's height above the 
water was the same for all the rudder angles. 

The forward flap was represented in Figure 11E. It 
started with a value of 0 degrees, reaching a peak of 13.9 
degrees to help the transient period. After 5 seconds it 
reached a steady state positive value of 10.8 degrees. This 
positive value means that the flap was oriented in the down 
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direction and therefore offering more resistance to the 
water, provoking more lift on the bow part of the craft. It 
can be concluded that the main roles of the forward flap are 
craft's lifting and to keep a stable pitch angle. 

In Figure 1 1F the stern mid foil segment was 
represented. After a transient period with a maximum value 
of -1.15 degrees, it reached a value of 2.62 degrees in 4.5 
seconds, as in the preceding case it was oriented in the 
down direction and again its main role is to help the lift 
motion of the craft. 

In Figure 1 1 G the stern starboard flap was 
represented. Its value was 0 degrees, before the rudder was 
introduced after 3 seconds. Then for a rudder angle of 30 
degrees, the flap had a transient period with a maximum 
value of -2.33 degrees, going to a steady state value of 
-1.86 degrees in 16 seconds. For a rudder angle of only 5 
degrees, the maximum negative value was only of -0.38 
degrees, going to a steady state value of -0.31 degrees in 
13 seconds. Having the starboard flap angles different 
values for different rudder angles, it was concluded that 
the stern starboard and port flaps have influence in the 
craft's turns, therefore a completely different conclusion 
in relation with the conclusion presented for the forward 
flap and mid aft flap. 

On the next groups of figures the same elements as 
before will be represented. Similar conclusions can be 
applied, therefore only the data obtained will be reported. 

2. Spe ed = 36 kno ts, rudder to starboard 

The pitch angle was represented in Figure 12A. It 
reached a steady state value of 0.01 degrees in 6.5 seconds. 

In Figure 12B were represented the roll angles for 
the different rudder angles, showing again the banking in 
the wrong direction. For a rudder angle of 30 degrees, the 
roll angle had a maximum negative value of -1.29 degrees, 
going to a steady state value of -1.12 degrees in 15 
seconds. For a rudder angle of 5 degrees, the maximum 
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negative roll angle was of -0.21 degrees, going to a steady 
state value of -0.18 degrees in 12.5 seconds. 

The heading angles represented in Figure 12C reached 
higher values due to a faster speed, increasing gradually 
with time, meaning that with the increase of speed the 
rudder effect on the craft's turning was larger. After 40 
seconds of run time with 30 degrees of rudder, the heading 
angle reached 8.5 degrees and for 5 degrees of rudder, the 
heading angle reached 1.42 degrees, clearly larger values 
than for the 30 knots speed. 

In Figure 12D the height above the water was 
represented, with a steady state value of 1.65 feet in 3.5 
seconds, after starting on 0 feet. 

The forward flap angle was represented in Figure 
12E. On the transient period it reached a maximum value of 
9.21 degrees, going to a steady state value of 7.7 degrees 
in 5.5 seconds. These values are smaller than for the speed 
of 30 knots, because the thrust is bigger and it is easier 
for the craft to lift. 

In Figure 12F the stern mid foil segment was 
represented. It reached a maximum value of -1.37 degrees in 
the transient period, jumping to a steady state value of 
1.74 degrees in 7.5 seconds. These values are also smaller 
than for the 30 knots speed by the same reasons of the 
preceding paragraph. 

In Figure 12G the different aft starboard flap 
angles were represented. For a rudder angle of 30 degrees 
the flap angle reached a maximum value of -1.94 degrees, 
going to a steady state value of -1.6 degrees in 16 seconds. 
For a rudder angle of 5 degrees, when the rudder was 
introduced the maximum value was -0.324 degrees, reaching a 
steady state value of -0.27 degrees in 13 seconds. 

3 • S peed =3 6 knot s^ rudder to port 

With the same speed but with the rudder oriented to 
the opposite board, port, there were no changes in the pitch 
mode, the only changes were verified with the roll angle, 
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angle, heading angle and aft starboard flap angle, 
represented in Figures 13A, 13B and 13C. The corresponding 
figures with the rudder to starboard were Figures 12B, 12C 
and 12G and comparing these two groups of figures, it can be 
observed that the roll angle, heading angle and aft 
starboard flap angle have exactly the same values but with 
opposite direction, because on the Y body axis direction the 
craft is perfectly sy mmetrical. 

4 . Sp ee d=40 knots 

In Figure 14A the several roll angles for different 
rudder angles were represented. When a rudder angle of 3C 
degrees was applied after 3 seconds the craft started 
banking in the wrong direction again, having a transient 
oscillation with a maximum value of -1.18 degrees. It went 
to a steady state value of -1.032 degrees in 15 seconds. 

In Figure 14B the different heading angles were 
represented. The craft went straight ahead and started 
turning when the rudder angle was introduced after 3 
seconds. The heading angles increased gradually with time, 
for a rudder angle of 30 degrees the heading angle reached 
the value 11.9 degrees after 40 seconds of run time. For a 
rudder angle of 5 degrees the maximum heading angle reached 
1.99 degrees in 40 seconds. 

In Figure 1 4C the height above the water was 
represented. It went to a steady state value of 0.628 feet 
above the water in 11 seconds. 

In Figure 14D the forward flap was represented. 
After a transient oscillation, with a maximum value of 6.9 
degrees, it dropped to a steady state value of 6.19 dgrees 
in 7 seconds. 

The aft starboard flap angles were represented in 
Figure 142. For 30 degrees rudder angle, the flap angle 
reached a maximum negative value of -1.78 degrees and went 
to a steady state value of -1.48 degrees in 16 seconds. For 
5 degrees rudder angle, the flap angle had a maximum 
negative value of -0.29 degrees and went to a steady state 
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value of -0.25 degrees in 11.5 seconds. 

5 • Sjree d=50 knots 

In Figure 15A the different roll angles were 
represented. For a 30 degrees rudder angle, the roll angle 
had a maximum value of -0.99 degrees and went to a steady 
state value of -0.88 degrees in 10.5 seconds. For a 5 degrees 
rudder angle, the roll angle reached a maximum value of 
-0.165 degrees and went to a steady state value of -0.147 
degrees in 8 seconds. 

The heading angles were represented in Figure 15B. 
For a 30 degrees rudder, the maximum heading angle was 18.3 
degrees in 40 seconds of run time. For a 5 degrees rudder, 
the maximum heading angle was 3 degrees after 40 seconds. 

In Figure 15C the forward flap angle was 
represented.lt went to a steady state value of 3.76 degrees 
in 3.5 seconds. 

In Figure 15D the different aft flap angles were 
represented. For a 30 degrees rudder angle, the starboard 
flap angle reached a maximum value of -1.51 degrees and went 
to a steady state value of -1.26 degrees in 14.5 seconds. 
For a 5 degrees rudder angle, the starboard flap angle 
reached a maximum value of -0.25 degrees and went to a 
steady state value of -0.19 degrees in 8 seconds. 
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Pitch angle in degrees 
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FIGURE 11 A 
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Boll angles in degrees 
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Reading angles in degrees 
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FIGURE 11C 
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Height above the water in feet 
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Forward flan angle in degrees 
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Stem start), flap angles in degrees 
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Roll angles in degrees 
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DISTURBANCE : RUDDJffi TO STARB. SPEED=36 KNOTS 
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Stem mid-flap angle in degrees 
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Stern starb. flap angle in degrees 
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Roll angles in degrees 
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Roll angles in degrees 
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Forward flap angle in degrees 
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Stem start), flap angles in degrees 
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Forward flap angles in degrees 
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B. EFFECTS OF A HELM STEP COMMAND 



Improving and completing the automatic control block of 
References 1 and 2 and based on the block diagram of Figure 
6.7 in Ref. 5, a simulated automatic control system was 
built, including the dynamic characteristics of the sensors, 
the dynamic characteristics of the compensation filters, 
inputs to the servo valves and the dynamic characteristics 
of the hydraulic actuators. 

There was an important limitation and perhaps a certain 
inaccuracy in the values used as corner frequencies of the 
compensation filters and the values used as gains for the 
control actuation servos, because References 1,2,3 and 4 
used the model PC(H)-’-1 for their studies and Ref. 5 used a 
more recent model, the PHM. Of course there was some lack, of 
information, but some interesting results and some data were 
obtained for future studies. 

As mentioned before, the main function of the forward 
flap was to lift the craft and keep it stable above the 
water in the pitch-heave mode. To confirm this, the helm 
doesn't actuate the forward flap, but only the starboard aft 
flap and the port aft flap, therefore these two last flaps 
will be responsible for the craft's turning in conjunction 
with the rudder and will also be responsible essentially for 
the craft's rolling and finally for the craft's pitching. 
The forward flap is only actuated according to its job by 
the forward accelerometer servo, by the height sensor servo 
and by the pitch angle servo. The two aft flaps are 
actuated by the respective accelerometer sensors, by the 
pitch angle sensor, by the roll angle sensor and by the 
manual helm control. The rudder servo is actuated by the 
roll angle servo, by a yaw rate gyro and by a heading 
command disengaged or engaged. In this last case it will 
stop and turn back the rudder when the desired heading is 
reached. 

Several figures will be shown with the roll angles. 
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heading angles, rudder angles, starboard and port flaps 
angles for different speeds and under the actuation of helm 
step commands of 150,120,90,60 and 30 degrees to the 
starboard side, at each speed. 

1 • Spe ed =30 knots 

In Figure 16A the roll angles corresponding to 
different helm angles were represented. The helm was moved 
to starboard after 12 seconds of straight ahead motion, then 
the craft started turning in the right direction, also to 
starboard. For a helm angle of 150 degrees the roll angle 
was 9 degrees after 24 seconds and for a helm angle of only 
30 degrees the roll angle was 3 degrees after 30 seconds. 

In Figure 16B the heading angles were represented. 
The craft started straight ahead and when a starboard helm 
was introduced, the craft started turning to starboard, 
reaching 21.2 degrees to starboard after 30 seconds for a 
helm angle of 150 degrees and only 3.4 degrees to starboard 
in the same amount of time for a helm angle of 30 degrees. 

Figure 16C shows the different rudder angles due to 
a helm step command. While the ship went straight ahead the 
rudder was also straight ahead, but when the helm was 
introduced after 12 seconds, the rudder turned to starboard 
on the craft's turning direction. Finally when the helm was 
removed after 28 seconds, the rudder went back to its normal 
position. For a helm angle of 150 degrees the rudder reached 
a value of 5.1 degrees to starboard and for a helm angle of 
30 degrees it only reached 0.78 degrees to starboard. 

In Figure 16D and 16E the port aft flap and the 
starboard aft flap were respectively represented. When a 
helm angle was introduced after 12 seconds, the craft 
started turning and banking. After those 12 seconds the port 
aft flap had positive values, (down direction) and the 
starboard flap had negative values, (up direction) . When the 
craft is turning to starboard, it must roll to starboard, 
therefore a lift on the port side must exist. The port flap 
being down, when the craft advances more force is exertad on 
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that flap by the water and there is a lift force on that 
side. The starboard flap is up and on that side there is a 
down force. When the helm goes to zero degrees after 28 
seconds, the craft rolls in the opposite direction and it 
can be verified by the figures that the flaps have also 
moved in opposite directions. 

2 • Speed=36 knots 

In Figure 17A the roll angles corresponding to 
different helm angles were represented. When the helm was 
introduced after 10 seconds, the craft started rolling to 
starboard. For a helm angle of 150 degrees the roll angle 
reached 6.2 degrees in 18 seconds and for a helm angle of 30 
degrees the roll angle reached 2.55 degrees after 27 seconds 
and it started banking in the opposite direction when the 
helm went again to 0 degrees. 

The heading angles were represented in Figure 173. 
When a helm angle of 150 degrees was introduced after 10 
seconds, the craft started turning to starboard, reaching a 
maximum heading angle value of 37.3 degrees after 28 
seconds and for a helm angle of 30 degrees, the craft 
reached a maximum heading angle value of 6.9 degrees after 
28 seconds. 

The rudder angles were represented in Figure 17C. 
For a helm angle of 150 degrees, the rudder angle reached a 
value of 5.37 degrees to starboard after 24 seconds, then it 
started going back in the opposite direction when the helm 
was removed. For a helm angle of 30 degrees the maximum 
rudder angle value was 1.44 degrees to starboard after 2*! 
seconds, then, as before, it started returning in the 
opposite direction. 

In Figure 17D and 17E the port and starboard flaps 
were represented. As before when the craft was banking to 
starboard, the port flap had positive values (down 
direction) and the starboard flap negative values (up 
direction) . Both flaps reversed the motions when the craft 
rolled back to its normal position after 24 seconds. 
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3 . 



S p ee <3=4 0 kno ts 

In Figure 18A the roll angles for different helm 
angles were represented. For a helm angle of 150 degrees, 
the roll angle reached a value of 9.38 degrees after 19 
seconds. For a helm angle of 30 degrees the roll angle 
reached a maximum value of 2.43 degrees after 23 seconds. 

In Figure 183 the different heading angles were 
represented. For a helm angle of 150 degrees the heading 
angle reached 26.3 degrees in 24 seconds and for a helm 
angle of 30 degrees the heading angle reached 12.12 degrees 
in 24 seconds. 

In Figure 18C the different rudder angles were 
represented. For a helm angle of 150 degrees the rudder 
angle had a maximum value of 5.83 degrees to starboard after 
21 seconds and for a helm angle of 30 degrees, the rudder 
angle had a maximum value of 1.16 degrees to starboard after 
21 seconds. 

In Figures 18D and 18E the port and starboard flaps 
were shown. Their movements were the same as before, 
according to the craft's rolling motion. 

4 . Speed=50 knots 

In Figure 19A the roll angles for different helm 
angles were represented. For a helm angle of 150 degrees 
the roll angle had a value of 9.45 degrees after 15 seconds 
and for a helm angle of 30 degrees the roll angle had a 
maximum value of 2.45 degrees after 20 seconds, then it 
started decreasing because the craft rolled in the opposite 
direction due to the removal of the helm command . 

In Figure 19B the heading angles were represented. 
For a helm angle of 150 degrees, the heading angle reached a 
value of 38.7 degrees to starboard after 25 seconds. For a 
helm angle of 30 degrees the heading angle reached a value 
of 11.1 degrees to starboard after 25 seconds. 

In Figure 19C the rudder angles were represented. For 
a helm angle of 150 degrees, the rudder angle reached a 
maximum value of 7.78 degrees after 21 seconds and then it 
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started decreasing, moving in the opposite direction due to 
the removal of the helm command. For a helm angle of 30 
degrees, the rudder angle reached a value of 1.54 degrees 
after 21 seconds and again it started decreasing. 

In Figures 19D and 19E, the port and starboard flaps 
were represented. Their motions are the same as before, but 
with different values. 
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Heading angles in degrees 
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Rudder angles in degrees 
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Roll angles in degrees 
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Rudder angles in degrees 
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Roll angles in degrees 
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C. RESULTS WITH THE SPLITTING OF THE FORWARD FOIL 

The forward foil was split in two parts to analyze how 
it improves or prejudices the craft's pitching, the craft's 
rolling and the craft's turning rate in the desired 
directions. 

Some changes in the computer simulaton program were 
necessary to achieve this trial. Thus, instead of 
considering a unique submergence for the forward foil, two 
different submergences. Seep and Secs, corresponding to the 
forward port foil and to the forward starboard foil were 
considered. These two foils were moved independently, thus, 
two different forward flap angles, Elvl and Elv3 were 
considered. Also two different forward foil drag 
coefficients and two different forward foil lift 
coefficients were introduced. Finally the block diagram of 
Ref. 5 was adapted to an extra forward foil, provoking a 
small change on the simulated automatic control system. 

With a helm step command to starboard, the craft 
started turning to starboard and it also started rolling to 
starboard. As was explained in part B of this chapter, the 
two forward flaps were moved in a manner similar to that of 
the aft port flap and the aft starboard flap, thus, the 
forward port flap went down provoking an upward lift force 
on the port side and the forward starboard flap went up 
provoking a downward force on the starboard side. Several 
runs with different speeds and different helm angles were 
done. In a general way, the roll angle in the desired 
direction was improved by 1055 and the turning rate was 
improved by 8%, but the pitch angle had values 3% smaller 
than the values obtained with a single forward foil, because 
with half of the forward foil in the upward direction, there 
are two lift forces, (one positive and the other negative) , 
prejudicing the craft's bow lift. 

Finally, the forward flaps were moved in opposite 
directions to those of the preceding paragraph, thus, the 
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forward port flap went up and the forward starboard flap 
went down. The results were worse than those of the preceding 
paragraph. 
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D. CALCULATION OF THE CORRECT ROLL ANGLE 



Many runs with the craft exposed to several 
disturbances were made. Different roll angles according to 
the existent conditions were obtained, but now a question 
arises : are these roll angles sufficiently near the optimum 
value ? 

A block called Calculation of the Correct Roll Angle 
was introduced in the computer simulation program to answer 
this question. Knowing the craft's turning radius for a 
certain speed and a certain helm angle, the craft's 
centrifugal force can be calculated. The craft's 
gravitational force in the Z direction was already known, 
then, the force perpendicular to the deck, called net force, 
can be calculated, by the gravitational and centrifugal 
forces. The angle between the gravitational force and the 
net force perpendicular to the deck is the correct roll 
angle. 

In Figures 20A and 20B the craft's trajectories in the 
X-Y plane for a speed of 40 knots and helm step commands of 
20 and 60 degrees respectively were represented. The craft's 
turning radius was calculated at the point where the heading 
angle was reversed 180 degrees in relation with the initial 
heading angle at the beginning of the run, since this will 
be approximately the steady state value. 

In Figure 20C the craft's turning path, for a helm 
command of 160 degrees, a speed of 50 knots (extreme 
conditions) and a run time of 156 seconds was represented. 
The turning path was approximately a circle, with an 
approximate turning radius of 410 feet. 

The values of the actual roll angle matched closely 
with the values of the correct roll angle, for the steady 
state period, although there were differences during the 
transient period. These values are represented in Table II. 

In Figure 20D the craft's turning path for a helm angle 
of 160 degrees and a speed of 30 knots was represented. The 



97 



turning radius was smaller, approximately 300 feet. 

The values of the actual roll angle matched again with 
the values of the correct roll angle, for the steady state 
period. These values are shown in Table III. 



TABLE II 



Roll angles for helm=160 degrees and speed=50 knots 



Time Actual roll angle 



(Seconds) 

17.3 
26.0 

34.7 

43.4 
52. 0 

60.7 

70.3 
78. 1 

86.8 

95.4 



(Degrees) 

13.98 

17.08 

18.33 

18.86 

19.14 

19.27 

19.36 

19.39 

19.39 

19.39 



Correct roll angle 
(Degrees) 

5.65 

10.22 

14.07 

18.19 

19.31 

19.39 

19.46 

19.52 

19.53 
19.53 
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TABLE III 



Roll angles for helm=160 degrees and speed=30 knots 



Time 

(Seconds) 
19. 2 

28.9 
38.5 
48.2 
57. 8 

66.9 
75. 2 



Actual roll angle 
(Degrees) 
10.77 
14.53 
14.88 
15.29 

15.47 

15.47 

15.47 



Correct roll angle 
(Degrees) 

7.32 
10. 5T 
13.84 
14.59 
15. 12 
15.72 
15.72 
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FIGURE 20A 
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TURNING PATH OF THE HYDROFOIL 
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TURNING PATH OF THE HYDROFOIL 
HELM=160 DEGREES SPEED=50 KNOTS 
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TURNING PATH OF THE HYDROFOIL 
HELM=l60 DEGREES SPEED=30 KNOTS 
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V* CONCLUSIONS AND RECOMMENDATIONS FOR FURTHER STUDIES 



The system was made stable in six degrees of freedom, 
with the following limitation : the dimensions, the foil 
areas, the weight and the thrust of the model PC(H)-1 were 
used in References 1 and 2 and in the computer simulation 
program. Due to a lack of information, because only very 
general information was available, certain values of the 
automatic control system of the model PHM, instead of the 
PC(H)-1, had to be used. 

A study of the craft turning characteristics was made 
in this thesis. The simulated automatic control system was 
improved and adapted to perform this study. 

Roll was controlled by differential displacements of 
the aft flaps, actuated by the helm control, while pitch was 
essentially controlled by the forward flap. 

The rudder has negative effects on the hydrofoil's 
rolling direction, which should be the same as the turning 
direction, due to the high speeds of operation. 

It was not worth while to split the forward foil. 

Recommendations for further studies include : 

1. Explore more deeply the advantages or the 
disadvantages of the splitting of the forward foil, 
combining it with the motions of the aft flaps and the 
rudder. 

2. Study the response of the craft in a sinusoidal sea 
and in a random sea. Design improved control systems for sea 
state operation. 

3. Simulate the transformation from hullborne operation 
to foilborne operation. 
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APPENDIX A 



SUMMARY OF COMPLETE EQUATIONS 

EQUATIONS FOR ACCELERATION IN BODY AXIS 

U = — F +RV-QW 
m X 

V » — F„+PW-RU 

m Y 

W = — F +QU-PV 
m Z 

» = J^t L-QR( I zz -I yy ) + (R+QP)I xz ) 

R = yi-[N-PQ(I yy -I xx )-(QR-P)I xz ] 

VELOCITY IN BODY AXES 

U = / Udt 

V = / Vdt 
W = / Wdt 
P = / Pdt 
Q *= / Qdt 
R = / Rdt 



(A-l) 

(A-2) 

(A-3) 

(A-4) 

(A-5) 



(A-6) 

(A-7) 

(A-8 ) 

(A-9) 

(A-10) 

(A-ll) 



TRANSFORMATION OF VELOCITIES FROM BODY TO EARTH .AXES 

U T? =UcosOcos^+V(cosi^sinO sin4)-sin^cos^ )+V/(cost|jsinOcos({) + sintj,'Sin<i) ) 

( A- 1 2 ) 

V^UcosOsini^+V (cos^cos4>+sini^sin0sin<|))+W(sin^sinOcos(?-cosi[/s.in<)) 
h ( A-13) 

W^=-UsinO+VcosOsin<^+VcosOcos(}) ( A-14 ) 



IO5 



<j> = P+^sin0 


( A-15) 


e = Qcos<}>-Rsin<J> 


(A-16) 


if> = (Qsin4>+Rcos<£ )cosQ + (4>-P)sin0 
POSITION IN EARTH AXES 


(A- 17) 


X E = / u E dt 


(A-18) 


y E = ! V E dt 


(A-19 ) 


z E = / w E dt 


(A-20) 


4> - f 4>dt 


(A-21 ) 


G = / Odt 


( A-22 ) 


= / tydt 

EXPANSION OF FORCE AND MOMENT EQUATIONS 


( A-23 ) 


F x - E f xif + E F xis + mg x +T x 


( A-24 ) 


Fy ^ Fy Jg + ITlgy 


( A-25) 


F z = z f zif + mg z 


( A-26 ) 



L ( F Z L Y^PF +( ' F Z L Y^SF ^ f y l z^ps ^ f y l z^ss ^ f y l z^cs 2/ ^ 



M E ( f z L X ) IF +E (F X L Z ) IF +E (F X L Z ) IS +T X L ZT 


( A-28 ) 


N = Z ( f y L X^IS 
EXPANSION OF GRAVITY TERMS 


( A-29) 


rn£ x = -mg sinG 


( A-30) 


mgy = mg cosGsincj) 


(A-31) 


mg^ = mg cos0sin<j> 


( A-32 ) 
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FOIL VELOCITY COMPONENTS IN BODY AXES 



Center Foil: 

u c - u+l zcf* q 

v c = v “ l zcf* p+l xcf* q 

w c ~ w " L xcp- ^ 



Port Foil: 



u p = u+l zpf‘ Q’W r 



v l zpf‘ p+l xpf* r 



Wp W+L YPF* P-L XPF* Q 



Starboard Foil: 



u s u+l zsf* q l ysf* r 



v s v l zsf* p+l xsf* R 



Wg W+L ySF *P L xsf *Q 



Mid Foil 



U M U+L XMF* Q 



V M V L ZMF’ P+L XMF* R 



V/ M W L XMF * Q 



(A-33) 
(A-34 ) 
(A-35) 



( A-36 ) 
(A-37 ) 
( A-3S ) 



( » nr\ \ 

/i-o: 7 ; 

(A-40 ) 
( A-41 ) 



( A-42 ) 
( A- 4 3 ) 
( A-44 ) 
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TRANSFORMATION OF WATER PARTICLE ORBITAL VELOCITY FROM EARTH 



TO BODY AXES 



U 



WI 



r WI 



W WI 



^EWI COS ^ COS ® +V EWI S i' n ^ COS ® - ^EWI S i. n €) (A-45) 

^EWI ^ cos ^ s i n 0 sin<J>-sim|;cos<J) )+V^ei ( simJ/sin0sin<t>+cos<{;cos<{> ) 

+ ^EWI COS ® s * n< ^ (A-46) 

^EWI ^ co s^ s in 0 cos^^sin^sincj) ) + ^£j ( sin^s in 0cos4>- cosies incf> ) 



+ W EW i cos 0 c ° s 4> 



(A-47 ) 



Each of these equations must be repeated for each foil. 



RELATIVE VELOCITY COMPONENTS 



U RI U I U WI 



V = V - V 
V RI I WI 



vi - w w 
W RI I WI 



One set of equations for each foil 



ANGLES OF ATTACK AND SIDE SLIP 



CXt = 



arctan W RI / U RI + “j.pjxED 
V T 



B t = arcsin 



RI 



/U 2 + V 2 +W 2 

/U RI RI RI 



( A-48 ) 
(A-49) 
(A-50 ) 



(A-51 )* 
( A-52 )* 



* One set of equations for each foil. Angles of attack and 
side slip are calculated in radians. 
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TOTAL RELATIVE VELOCITY AT A PARTICULAR FOIL OR STRUT 



V - /u’ + vj+w‘ 


( A-53 ) 


HYDRODYNAMIC FORCES IN WATER AXES 





(Each foil equation adn strut equation must be repeated for 
foil and strut, respectively) 

Foil Lift: 



F L = 4pv ! a p c f 
Foil Drag: 


(A-54) 


F DF = ^PV 2 A p C p 
Strut Drag: 


(A-55) 


F DS - if>V ’ A S C DS 
Strut Side Force: 


( A-56 ) 


F s = ipV 2 A s C s 


(A-57 ) 


TRANSFORMATION OF FOIL FORCES FROM WATER AXES TO BODY 


AXES 


F xif = -F IDF cosacosB+F LIF sina 


( A-58 ) 


F YIF = “ F DIF Slne 


(A-59) 


F ZIF " -F DIF sinacosB-F LIF cosa 


( A-60 ) 
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// EXEC DSL , PARM.C='NOSOURCE ,NOMAP* 

//DSL. INPUT DD * 

D COMMON M ATR I X ( 53 9 ) 

CONTRL FINTIM=30.0,DELT=0.025,DELS=0. 1 
INTEG RKSFX 

INTGER I FLAG, I FL AG1 , NP L OT , NUM , NCUR 
CONST N FLOT= 1 ,NCUR= 1 , IFLAG=0 

CONST WE I GHT = 262 52 7, IXX=60I200, I YY=4286000, . . . 

I Z Z = 4l 74000, G = 32. 2 

CONST ACF=32. 7 ,APF=36.3 , ASF=36.8, AMF=76.0, . . . 

PI=3. 141593 

CONST LXCF =41 . 5 9 , LXCS = 4 1 . 5 9 , LX PF = - 1 6 . 9 8 , LXPS = - 16 . 98 , .. . 
LXSF=-16.93,LXSS=-16. 9 8 , LXMF =-l 6 . 98 , L YCP=-5 . 0 , LYCS=5.0, . . . 
LYPF=-12.42 ,LYPS=-8.G, LYSF= 12.42, LYSS=8.0,LZCF= 15.5, . . . 

LZCS = 13.25,LZFF = 17.14,LZPS=14.0,LZSF=17.14,LZSS=14.0,. . . 

L ZM F= 1 7 .14,LXT=-16.9, LZT=17. 7 » LXH=68 . 0 
CONST K1=5.0,K2 =-9 8.0,K3=-82.0 

CONST K4=l . 5 ,K5 = 10.0, K 6=3.5 , K7= 1 . 42 , K8= 1 .5 , K9 =- 2 . 0 , . . . 

K10=0.44,Kll=-13.0,Kl2=3f0,K13=3.0,K14=0.014 

CONST R ST = 0 . 0 

STORAG ICK3) ,NU1( 3) , DEN 1 1 4 ) , IC2( 2) » NU2 ( 2 ) ,DEN2(3), ... 

I C3 ( 3 ) , NU3 ( 2 ) , DEN3 ( 4 ) 

TABLE IC1(1-3)=3*0.0,NU1 ( 1-3 } = 1 .0 , 3. 0 1 , 0 . 03 , . . . 

DEN l ( 1-4)= 1.0, 110. 1, 3011.0,3 00. 0,IC2( 1-21*2*0.0,. .. 

NU2( 1-2)=1.0,1.0,DEN2( 1 -3 ) = 1 . 0 , 87 . 7 , 1 885 .0 , IC 3 ( 1- 3 ) =3*0 . 0, . . 

NU3( 1-2) =30. 0,5 8.5, DEN3( 1-4) =1.0, 30. 25, 0.0 15, 0.45 

PARAM UZER0=40 . 0 

PARAM RUDR=0 . 0 

PARAM HELMR=2 .79 

PREPAR .1,RTIME»PHI1,PSI1,XE,YE 

GRAPH XE , YE 

PRPLCT 

LABEL RESPONSES WITH HELM STEP COMMAND 

PRI NT 0. 5, R TIME, PH I , PH 1 3, PH II , PH 1 4 , PS I , PS I 1,XE ,YE ,UE,VE , . . . 
RGIR,FGRAZ,FCENT,FNET 

INI TIAL 

* 

CENSE = 2 .0 

U I NIT=U ZERO* 2000. 0*3.0*1.0/3600.0 

U=U I NIT 

TX = 3960 .0 

ALM=0. 074269 

ALS=0. 074269 

ALP=0. 074269 

ALF=0. 037111 

PREF=0 .0 

TREF=0 . 0 

THETA2=0.0 

PH 12 = 0 . 0 

PS I 2=0 . 0 

MASS= WE I GHT/G 

CCNV=P 1/180.0 

HREF= 1 1 .5 

HR EF2 = - 1 1 .5 

LENGTH= 115.75 

FLAG=0. 0 

I FL AG I = 0 

E PS=2 5 . 0 

I F ( I FLAG . EQ. 0 ) CALL S ETUP ( F LAG ) 

I FLAG = 1 

IF(FLAG.EQ.l) CALL ENDRUN 

DERIVATIVE 

* 

NOSORT 

* ANGLES OF ATTACK AND SIDE SLIP 

* 

VRC=V-LZCF*P+LXCF*R 
VRS=V-LZSF*P+LXSF*R 
VRM=V— LZMF*P + LXMF*R 



no 






VRP=V-LZPF*P+LXPF*R 

WRP=W+LYPF*P-LXPF*Q 

WRS=W+LYSF*P-LXSF*Q 

WFM=W-LXMF*Q 

WRC=W-LXCF*Q 

ALFC=WRC/U+ALF 

ALFA P=W RP/U+ALP 

ALFAS=WRS/U+ALS 

ALFAM=WRM/ IM-ALM 

BETAC=VRC/U 

BETAP=VRP/U 

8ETAS=VRS/U 

CALCULATION OF EULER ANGLES 

PHID1=P 
TFETD1 =Q 
PSID1 =R 



THETA=I NTGPL ( THE T A2 , T HETD 1 J 
PHI=INTGRL(PHI2,PHIQ1) 

PS 1 = 1 NTGRL ( PS 12, PS ID1 ) 

WE=-U*S I Nl THETA) *V*COS ( THETA)*S IN < PHI ) *W*COS< THETA ) . .. 
*COS< PHI ) 

ZE=INTGRLIHREF2,WE) 

ZE1=ZE-HREF2 



SUBMERGENCE TERMS 

SECF=ZE-LXCF*THETA*LZCF*COS(PHI > 
SEPF=ZE-LXPF*THETA+LYPF*SIN(PHI ) +LZPF*COS < PHI ) 
SESF=ZE~LXSF*THETA+LYSF*SIN( PHI ) *LZSF*COS ( PHI ) 
SEMF=ZE-LXMF*THETA + LZMF*COS(PHI ) 

HYDRODYNAMIC COEFFICIENTS 

HELMR1=STEP(RST) 

HELM=HELMR*HELMRl 
RUD1 1 = K7*PHH-K9*R 
RUD=REALPL(O.Q, 0.159, RUD11) 

ELV1=EL V*CCNV 
ELVD=LIMIT(-5.0,5.0,ELV1) 

FLA D= FLAP* CON V 

FDEG=L1 MIT(-5.0,5.0,FLAD) 

ARPD=ARP*CONV 

ARP 1 = LI MIT (-5 .0,5 .0, ARPD) 

ARSO=ARS*CONV 

ARSli=LIMI T< - 5.0 ,5 .0, APSO) 

CALL INTERPt ALF C , E LV D, CGEFC , 2 ) 

CALL INTERPI ALFAP , ARP1 ,C0EFPX»3) 

CALL INTERPI ALFAS, ARS11 ,COEFSX t 3) 

CALL INTERPI AL F AM , FLA P , COEF MX , 1 ) 

CALL INTRP1(BETAC,SECF,CFDS»2,THETA,THETA1 ) 
CALL INTRPlt BET AP , SEPF ,CPDS, I , TH ET A , THE T A I ) 
CALL INTRPlt BETAS, SESF.CSDS, 1, THETA, THETA I ) 

FORCES IN THE WATER AXES 

FGRAX=-MASS*G*S IN (THETA) 

FGRAY=MASS*G*SIN ( PHI ) 

FGRAZ=MASS*G*COS< PHI 1 

CCS =-0.043 75* BETAC-t-0 .04375*RUD 

ACCS=CCS*4.25*SECF 

FSC=0. 5 *DE NS E*i U-*2) * ACC S 

FDSC=0. 5*DENSE*( U**2) *CFDS 

PCS=- 2 . 1 5* 3ET A P 

APCS=PCS*4.25*SEPF 

FPS=0. 5*DENSE*l'J**2 )*APCS 

FD PS= 0 . 5* DENSE*! U**2) *CPDS 
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SCS=-2. 15*BETAS 
ASCS=SCS*4.25*SESF 
FSS=0.5*DENSE*( U**2)*ASCS 
FCSS=0. 5* DENSE*! U**2) *CSDS 
CLCF= 5 . l5*ALFC+0 .036*ELV+0. 13 
FLC = 0. 5*DEiNSE*(U**2) *ACF*CLCF 
FDC=0.5*DENSE*( U**2 ) *ACF*COEFC 
CLPF=4. 10*ALFAP<-0.0L8*ARP 
FLP=0. 5*DENSE*<U**2 >*APF*CLPF 
FCP=0.5*DENSE*<U**2) *apf*coefpx 
CLSF= 4. 10*ALFAS+0.018#ARS 
FLS=0.5*DENSE*(U**2)*ASF*CL$F 
FDS=0.5*DENSE*( U**2) *ASF*COEFSX 
C LMF= 5 . 15*ALFAM+0.036*FLAP 
FLM=0.5*CE.NSE*( U**2 )*AMF*CLMF 
FDM=0.5*DENSE*( U**2)*AMF*C0EFMX 
F L C Z = 0 . 5*DENSE*( U**2J *ACF*CLCF 
FCCZ=0. 5*DENSE*(U**2> *ACF*COEFC 
FLPZ = 0. 5* DEM Sc* ( U**2) *APF*CLPF 
FCPZ=0.5*CENSE*(U**2 ) *APF*CCEFPX 
FLSZ=0. 5*DENSE*( U**2) *ASF*CL$F 
FDSZ=0 . 5*DENSE*( U*^2) *ASF*COEFSX 
FLMZ=0.5*DENSE*(U**2)*AMF*CLMF 
FDMZ=0. 5*DENSE*( U**2I *AMF*COEFMX 



TRANSFORMATION TO BODY AXES 

FXCF=-FCOALFC*FLC 

FXPF=-FDP+ALFAP*FLP 

FXSF=-FDS+ALFAS*FLS 

F X M F= -F DM <■ AL F A M* F L M 

FXCS=-F DSC-BET AC*FSC 

FXPS=-FO PS-BET AP*F PS 

FXSS=-FCSS-BETAS*FSS 

FYCS=-BETAC*FDSC*FSC 

FYPS=-BETAP*FDPS+FPS 

FYSS=-BETAS*FDSS<-FSS 

FZCF=-FDC#ALFC-FLCZ 

FZPF=-F DP Z*ALF A P-FLPZ 

FZSF=-FCSZ*ALFAS-FLSZ 

F ZMF=-F OMZ-Y-ALFAM-FLM Z 

LZCS=LZCF-0.4*SECF 

LZPS=LZPF-0.4*SEPF 

LZSS=LZSF-0.4*SESF 

FXCNE=FXCF+FXPF+FXSFfFXMF+FXCS*FXPS+FXSS 

* 

SUMMATION OF FORCES AND MOMENTS 

FX=FXONE + FGRAX+-TX 

UDOT= ( 1 . / M A S S ) * F X - Q* W 

IF< IFLAG1 .EQ.O) UDOT^O.O 

FY=FYCS+FYPS+FYSS+FGRAY 

FZ=FZCF*FZPF*FZSF+FZMF+FGRAZ 

L=FZ PF*LY PF*FZSF*LYSF-FYPS*LZPS-FYSS*LZSS-FYCS#LZCS 
MFWD=FZCF*LXCF 

K=-MFWD-FZPF*LXPF-FZSF*LXSF-FZMF*LXMFf. . . 

FXCF*LZCF + FXPF*LZPF+FXSF*LZSF + FXMF*LZMF«-FXCS.. . 

*LZCS+FXPS*LZPS+FXSS*LZSS+TX*LZT 

N=FYCS*LXCS+FYPS*LXPS+FYSS*LXSS 

VELOCITY INTEGRALS 



SORT 

VDCT= ( 1 . /MASS ) *FY-R*U 
WDCT = ( 1 ./MASS)*FZ + Q*U 
PDCT= ( 1 .0/ IXX)*L 
CDOT= ( 1 .0/IYY)*M 
RDCT= ( 1 .O/IZZ) *N 
U = I NTGR L ( U I N IT » UDOT ) 
V=INTGRL(O.G»VDOT ) 
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to=INTGRL(O.Q , WDOT ) 

P=INTGRL(0.0,P00T) 

Q=INTGRL(O.0,QDOT) 

R=INTGRL(0.0,RD0T) 

* SENSED ACCELERATIONS 

❖ 

AZC=FZ/MASS-LXCF*QDOT 

AZS=FZ/MASS+LXSF*QDOT+LYSF*PDOT 

AZP=FZ/MASS+LXPF*QDOT-LYPF*POOT 

* AUTOMATIC CONTROL BLOCK 

* * 
HACT=-ZE+LXH*THETA 

DELH=HREF-HACT 

DELHC = REALPL( 0. 0,0.04, DELH ) 

DETH=TR EF-THETA 
FLAP=K2*DETH 
B=AZC<-K14*{ PHI**2 ) 

AZCC=TRNFR(2,3,IC3,NU3,DEN3,8) 

A ZSC=RE ALP L( 0.0,0.033 , AZS) 

AZPC=REALPL( 0. C, 0.033 , AZP) 

THETAC=TRNFR (2,3 ,IC1 ,MU1 , DENI , 1 000 . 0*TH E T A ) 

PHIC=TRNFR ( I ,2, IC2,NU2,DEN2, 1865 . 0*PHI I 
HELMC = R EALPL ( 0 . 0 , 1.0, HELM) 

ELV11=K1 l*THETAC-«-K12*THETA<-K13*DELHC+K4*AZCC 
£LV=REALPL(0.0,0.I59,ELV11 ) 
ARS2=K4*AZSC+K5*THETAC*K6*PHIC-K10*HELMC 
ARP2=K4*AZPC+K5vTHETAC-K6*PHIC*K10*HELHC 
AFS= REAL PL (0.0,0. 159, ARS2) 

ARP=R EALPL (0.0, 0.15 9, ARP2) 

TFETA3=( THETA-THETA2) *130. 0/P I 
PHI1=( PHI- PH 12)^180. O/PI 
P SI I = < PS I -PS 12) - 1 80.0 /PI 
RUD1=RUD*180.0/PI 

❖ 

* CALCULATION OF THE CORRECT ROLL ANGLE * 

❖ 

UE=U*CCS( THETA )*C OS ( PS I ) +V*(COS( PSI) *S IN ( THETA) *.. . 
SIN( PHI) -SIN( PS I ) *COS (PHI ) ) <-W* ( C OS ( P S I) * S I N ( THET A ) . . . 
*CGS( PH I H-SI N( PS I )*SIN( PHI ) ) 

VE=U*COS( THETA) *S IN( PSI ) f V* ( COS ( PSIJ *COS ( PH I ) +• . • 

SIN( PS I )*$IN( THETA ) *S INC PHI ) ) +W*( SIN(PSI )*S IN (THETA) .. 
*CCS( PHI ) -COS ( PSI)*SIN(PHI)) 

YE=INTGPL(0.1,UE) 

XE=INTGRL<0.1 ,VE) 

RG I R= ( ( XE**2)<-( YE**2) ) /XE 
FCENT=MASS*< ( ( UE** 2 ) + ( VE**2 ) ) /RGIR) 

FNET=SQRT( (FGRAZ**2 ) M FCENT**2) ) 

PHI 3 = AT A N2 ( EC ENT , FGRAZ i 
PH 14= PH 13*180 . 0/P I 



* CALCULATION OF REAL TIME 

❖ 



F T I M E= ( LENGTH/U INIT)*TIME 

* 



DYNAMIC 

* 



SAMPLE 

CALL DRWG ( 1 , 1 , X E , YE ) 

TERMINAL 

CALL END RW ( N PLOT ) 

END 

STOP 

FORTRAN 

C 

C 

C SUBROUTINE SETUP 

C THIS SUBROUTINE FORMS THE DATA POINTS COMPRISING 

C THE CURVES OF FOIL AND STRUT DRAG COEFFICIENTS 

C AND SIDE SLIP COEFFICIENTS INTO MATRACc S . 

C 
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SUBROUTINE SETUP(FLAG) 

COMMON COEFMX,CO£FCX,COEFOX,CFDSl ,CADS1 
DIMENSION COEFCX (3*33) ,COEFOX( 3 , 33} » COE F MX ( 3 , 33 I t 

1CFDSI ( II , II ) , CAD SI ( 1 L ,11 ) , INFO( 8) 

REAL INFO 
IC = 0 
I FLAG = 0 

I FLAG 1 = 0 
M AT RI X= 1 
J = 0 

K = 1 

I CCUNT=0 

2 R E AD( 5 » 100) INFO 
100 FORMAT ( 7 FI 0. 8 ) 

I I = 1 

I F ( IFLAG1.E0.1) GO TO 51 
GO TO (4,5,6,7,81, MATRIX 

4 CONTINUE 
DO 1 1=1,8 
J = J + 1 

IFU.GT.3) K=K + 1 
IFU.GT.3) J = 1 
CCEFOXl J,K)=INFO( I ) 

ICCUNT= ICOUNT+1 
IC=IC+1 

I F ( I COUNT . EO .99 ) GO TO 3 
1 CONTINUE 
GO TO 2 

3 11 = 1 

IFU1.LE.8) 11 = 11+1 
J = 0 
K = 1 

I CCUNT=0 
MATRIX=2 

5 IF( II .GT.8) GO TO 2 
DO 9 1=11,8 

J = J + 1 

IFU.GT.3) K=K + 1 
IFIJ.GT.3) J=1 
COEFCXt J,K) = INFQ{ I ) 

I COUNT= IC0UNT+1 
IC=IC+1 

IF ( I COUNT « EQ.99 ) GO TO 10 
9 CONTINUE 
GC TO 2 
10 11 = 1 

IFU1.LE.8) 11 = 11+1 
J = 0 
K= I 

ICCUNT=0 
MATRI X= 3 

6 I F ( II .GT.8 ) GO TO 2 
DO 11 1=11,8 

J = J+1 

I F ( J . GT . 3 ) K=K+ 1 
I F < J . GT .3 ) J = 1 
CCEFMX ( J , K )= INFO ( I ) 

ICCUNT= ICOUNT+1 
IC=IC+1 

IF < ICOUNT. EQ.99) GO TO 12 

11 CONTINUE 
GO TO 2 

12 11=1 

IFU1.LE.8) 11=11 + 1 
J = 0 
K= 1 

I COUNT = 0 
MATRI X=4 

7 I F ( I 1 .GT.8) GO TO 2 
DC 13 1=11,3 

J = J + 1 
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IF(J.GT.ll) K = K + 1 
IFU.GT.ll) J=1 
CADS1( J » K ) = J NFO ( I ) 

I C CUNT = ICGUNT + 1 
IC=IC+1 

I F ( I COUNT . EG . 121 ) GO TO 15 

13 CONTINUE 
GC TO 2 

15 11=1 

I F ( II . LE • 8 ) 11 = 11+1 
J = 0 
K= 1 

I CCUNT = C 
M ATR I X= 5 

8 I F ( 11 .GT.8 J GO TO 2 
DC 14 1=11,8 
J = J + 1 

IF(J.GT.ll) K=K + 1 
IF(J.GT.ll) J=1 
CFDS1 ( J , K) = 1 NFO( I ) 

I COUNT = IC GUNT + 1 
IC=IC+1 

IF( ICGUNT.EQ.121 ) GC TO 16 

14 CONTINUE 
GO TO 2 

16 I FLAG 1 = 1 

50 I F LAG= 1 

I F ( ( IFLAG.EQ.l ) . AND. ( I FLAG1 .EG. 1 ) ) GO TO 52 

51 WRITE(6,104) ICOUNT, MATRIX, IC 

104 FORMAT ( ’ ‘,’JOB F A I LURE • , L X , • C OUNT= * , I 6 , IX , • MAT RJ. X= * 
116, IX, ' NO. PTS= 1 , 16) 

FL AG= 1 
RETURN 

52 WRITE ( 6, 105 ) IC 

105 FORMAT ( * * , • JOB COMPLETE 1 X, • POINTS ST0RED=‘,I6) 

RETURN 

END 

FORTRAN 



SUBROUTINE INTERP 






SUBROUTINE INTERP(ALFAC,FLAP»COEF,M) 

COMMON COEFMX, CO EFCX , C OEFOX 

DIMENSION COEFMX (3,33) , COEFCXl 3 , 33 ) , COE FOX ( 3, 33 ) 

DATA CONDEG/57. 295773/ 

ALFA=ALFAC*COND£G 

I F ( ALFA . GT. 10.0) WRITE(6,100)ALFAC,M 
1F( ALFA. GT. 10.0) ALFA=10.0 

100 FORMAT ( *0‘ , ' ERROR. ALFA EXCEEDED L I M ITS . A LF A=» , E 1 4. 8 , 
12X, • M=' ,11) 

I F( ALFA .LT.-6.0J WRITE(6,100) ALFAC, M 

IF(ALFA.LT .-6.0) ALFA=-6.0 

AL PHA = 2 . 0* ALF A 

IF( ALPHA. LT.0.0 ) GO TO 1 

IF< ALPHA. GT. 0.0) GO TO 2 

I A=1 3 

WT A=0. 0 

GC TO 3 

1 I A 1=- AL P HA 
I A = 1 3- I A 1 
ALPHA1 = - I A 1 
WTA=ALPHA- ALPHA1 
GO TO 3 

2 I Al= ALPHA 
I A=1 3+ I A1 
AL PHA 1= I A1 
WTA=ALPHA-ALPHA1 

3 IF( FLAP. LT. -5.0) WRITE(6,101) FLAP 
IF(FLAP.LT.-5.0) FLAP=-5.0 
IF(FLAP .GT.5.0) WRITE(6,101) FLAP 
IF(FLAP.GT .5.0) FLAP=5.0 
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101 FORMAT* •0* TERROR. FLAP EXCEEDED LI MI TS . FLAP = ' 
IF(FLAP)4,5,6 
A I B= 1 

WTB=-( ( 5.0+FLAP) /5.0) 

GO TO 7 

5 18 = 2 
WTB=0.0 
GO TO 7 

6 I B = 3 

WTB={ 5 .0-FLAP ) / 5 .0 

7 I F(WTA) 8,9,10 

8 I A2= I A- 1 
GO TO 11 

9 I A 2= I A 
GC TO II 

10 I A2=I A+ 1 

11 IF(WTB) 12,13,12 

12 I B2= 2 

GG TO 14 

13 IB 2= IB 

14 GC T0(15,16,17) ,M 

15 C0EF1 = C0EFMX ( 13, I A > 

C0EF2=C0E FMX ( 182, I A) 

GC TO 18 

16 CCEF1=CCEFCX( IB, IA) 

CCEF2 = C0EFCXII B2 , IA) 

GC TO 18 

17 C0EF1=C0EF0X( IB, IA ) 

C0EF2=CCEFQX ( IB2, IA) 

18 C0EF3=( C0EF2-C0EF1 ) * ABS( WTB ) +COEF 1 
GO TO ( 19,20,21) ,M 

19 CCEF1P = C0EFMX( IB, IA2) 

C0EF2 P=COE FMX ( 162, IA2) 

GC TO 22 

20 C0EF1 P=COEFCX( IB, I A2 ) 

CC£F2P=C0EFCX( IB2, IA2) 

GC TO 22 

21 CGEF1 P=COEFOX( IB, I A2) 

COEF ; ZP = COEFOX( IB2, 1A2) 

2 2 CCEF3P= (CCEF2P-C0EF1P ) *ABS( WTB ) 4-CCEF1P 
CCEF= (COEF3P-COEF3 ) * AB S ( WTA ) <-C0EF3 
RETURN 
END 

FORTRAN 



SUBROUTINE INTRP1 



E14. 8) 



SUBRGUT INF. INTRP1 ( BET A1 , SOB ,C OEF S , M, THE T A , THETA1 ) 
COMMON JUNK, CFDS 1 , CADS 1 

DIMENSION JUNK (297) , CFDS 1 ( 1 1 , 1 1 ) ,CADS1( 11,11) 

DATA CONDEG/57. 295773/ 

SUB=SOB 

BET A= BETA l^CONDEG 

IF(BETA.GT.15«0) KRITE(o,100) 8 E T A 1 , M 
I F ( BE TA «> GT .15.0) BETA=15.0 

100 FORMAT ( * 0 * , * ERROR . BET A EXCEEDS LIMITS. BETA=* »E14.3,2X, 
1' M=' , II ) 

I F ( BE TA • LT 1 5 . 0 ) W P.I TE (6 , 1 00 > BETA1 ,M 
IF(BETA.LT.-15.0 ) BETA = -15.0 
IF(BETA.LT.O.O) GC TO 1 
IF(BETA.GT.O.O) GO TO 2 
I A = 6 
W T A = 0 . 0 
GC TO 3 

1 I A 1 = -BE TA 
I A = 6— i A 1 
BET AG=- I A 1 
WTA=BETA-BETAG 
GC TO 3 

2 I A 1 = 8 E T A 
I A=6+ I A 1 
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BETAG=I A1 
WTA=BETA-BETAG 

3 IF(SUB.LT.O.O) SUB=0.0 
IF(SUB.GT.IO.O) SUB= 10 . 0 
IF(SUB.GT.O.O) GO TO 5 
I S= 1 
WT S = 0 . 0 
GO TO 7 

5 I S1=SUB 
I S= I S 1 + 1 

see 1 = i s i 

WT S=SU8-SUB1 

7 I F( WTA) 8,9, 10 

8 I A2= I A- 1 
GC TO 11 

9 I A2 = I A 
GC TO 11 

10 IA2=XA+1 

11 IF(WTS.GT.G.O) GO TO 12 
I S2= I S 

GC TO 15 

12 I S2= I S+ 1 

15 I F ( M » E G . 2 ) GO TO 16 
C0EF1=CADS1 ( IS, I A > 

CCEF2=C ADS11 2S2 , I A) 

GO TO 18 

16 CCEF1=CFDS1( IS, I A I 
CCEF2=CFDS1( IS2, IA) 

18 C0EF3=( C0EF2-C0EF1 ) *WTS+C0EF1 
I F ( M • EQ . 2 ) GO TO 20 

19 CCEFi P=CA0S1 ( IS, I A2) 

CCEF2P = CAD$1(IS2, IA2) 

GC TO 22 

20 COEF 1P=CFDS1 ( IS, IA2) 

CCEF2P = CF0S1 1 IS 2 , IA2) 

22 CCEF3 P= (C0EF2P-C0EF1P)*WTS*C0EF1P 
COEFS= (C0EF3P-C0EF3)*ABS< WTAJ +C0EF3 
THETA1=THETA*C0NDEG 
RETURN 
END 

//L.SYSUT1 DD SPACE=< 3500, ( 100, 10) ) 
//L.SYSLMOD DO SPACE= ( 3500, ( 100, 10, 1 )) 

/ / G . FT 06F001 DD SPACE = ( C YL , < 31) 
//G.AUXCATA DD * 

0.018 
0.024 
0.019 
0.017 
0.024 
0.0135 
0.017 
0.023 
0.016 
0.016 
0.022 
0. 0165 
0.016 
0.022 
0.017 
0.015 
0.021 
0.017 
0.015 
0.020 
0.017 
0.014 
0.020 
0. 0175 
0.014 
0.019 
0.018 
0.013 
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o. oia 
0.0 18 
0.013 

o.oie 

0.019 
0.013 
0.018 
0.020 
0.013 
0.013 
0.020 
0.013 
0.019 
0.021 
0.013 
0.020 
0.022 
0.0135 
0.020 
0.023 
0.014 
0.022 
0.025 
0.014 
0.023 
0.026 
0.015 
0.024 
0.028 
0.016 
0.026 
0.030 
0.0175 
0. 0275 
0 . 0^2 
0.020 
0 . 03 C 
0.035 
0.021 
0.032 
0.0375 
0.023 
0.034 
0 . 040 
0 . 025 
0.037 
0.044 
0.028 
0.040 
0.0475 
0.030 
0.044 
0.052 
0. 034 
0.047 
0.055 
0.036 
0.051 
0.059 
0.040 
0.054 
0.068 
0.044 
0.058 
0.067 
0.047 
0.063 
0.070 
0.050 
0.065 
0.074 
0.025 
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0.0 26 
0.027 
0.024 
0.025 
0.026 
0.023 
0.023 
0.0 24 
0.020 
0.021 
0.022 
0.018 
0.019 
0.021 
0.0165 
0.0185 
0.020 
0.015 
0.017 
0.013 
0.014 
0.015 
0.0175 
0.013 
0.0145 
0.013 
0.013 
0.0 14 
0.017 
0.011 
0.013 
0.018 
0.010 
0.012 
o.oie 
0.008 
0.012 
0. 019 
0.007 
0.012 
0.020 
0.006 
0.012 
0.021 
0.005 
0.013 
0.023 
0.005 
0.0135 
0.028 
0.005 
0.014 
0.026 
0.0055 
0.015 
0.027 
0.006 
0.016 
0 .G 30 
0.0065 
0.0175 
0.0325 
0.007 
0.020 
0.035 
0.008 
0.022 
0.038 
0.010 
0.0245 
0.043 
0.012 
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0.0265 
0.047 
0.014 
0.0245 
0.054 
0.016 
0.033 
0.060 
0.018 
0.037 
0 . 07 C 
0.021 
0,042 
0.078 
0.028 
0.048 
0.084 
0.026 
0.056 
0.093 
0.032 
0.067 
0 . 10 
0.036 
0.08 
0 . 1 1 
0.021 
0.021 
0.0215 
0 . 02 C 5 
0.0205 
0.021 
0.020 
0.020 
0.0205 
0.018 
0.018 
0.019 
0.017 
0.0175 
0.018 
0.0165 
0.017 
0.0175 
0.015 
0.0155 
0.016 
0.014 
0.0145 
0.0155 
0.0135 
0.014 
0.015 
0.0125 
0.013 
0.0145 
0.01C 
0.012 
0.015 
0.0095 
0.01 1 
0.015 
0.008 
0.0105 
0.016 
0.0075 
0. 01C 
0.017 
0.006 
0. 0105 
0.0195 
0.005 



0.011 
0.022 
0.0045 
0 . 0 L 2 
0.024 
0.0045 
0.013 
0.028 
0.005 
0.0145 
0.031 
0.005 
0.0165 
0.035 
0.006 
0.020 
0.041 
0.008 
0.023 
0.048 
0.010 
0.027 
0.055 
0.012 
0.031 
0 • 064 
0.015 
0.036 
0.078 
0.019 
0.044 
0.093 
0.022 
0.051 
0 . 10 
0.028 
0.060 
0 . 10 
0.035 
0.075 
0 . 10 
0.039 
0.084 
0.10 
0.045 
0.092 
0 . 10 
0.05 
0.095 
0.10 
0.056 
0.10 
0.10 
0.25 
0 . 30 
0.345 
0.41 
0.46 
0.52 
0.56 
0.62 
0.65 
0.70 
0.78 
0.25 
0 . 295 
0.34 
0.38 
0.425 
0.47 
0 . 52 
0.55 
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0.60 
0.63 
0.68 
0 . 25 
0.29 
0.325 
0.36 
0.40 
0.43 
0.47 
0.50 
0.545 
0.575 
0.62 
0.25 
0.28 
0.32 
0.345 
0.3 7 
0.41 
0.44 
0.47 
0.50 
0.54 
0.57 
0.25 
0.28 
0.31 
0.34 
0 . 36 
0.40 
0.43 
0.45 
0.48 
0.51 
0.54 
0.25 
0.28 
0.31 
0 . 34 
0.36 
0.39 
0.4 2 
0.45 
0.47 
0.50 
0.53 
0.25 
0.28 
0.31 
0.34 
0.36 
0.40 
0.43 
0.45 
0.48 
0.51 
0.54 
0.25 
0.28 
0.32 
0.345 
0.37 
0.41 
0.44 
0.47 
0.50 
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0.325 



122 



0.36 
0.40 
0.43 
0.47 
0.50 
0 . 545 
0.575 
0.62 
0.25 
0 . 295 
0.34 
0.38 
0.425 
0.47 
0.52 
0.55 
0.60 
0.63 
0.68 
0.25 
0.30 
0.345 
0.41 
0.46 
0.52 
0 . 56 
0.62 
0.65 
0.70 
0.78 
0.05 
0.11 
0.17 
0.23 
0.29 
0.35 
0.40 
0.46 
0.52 
0.58 
0 . 64 
0.05 
0.10 
0.15 
0.195 
0 . 245 
0.29 
0.34 
0.38 
0.44 
0.48 
0.54 
0.05 
0.095 
0.14 
0.17 
0.21 
0 . 25 
0.28 
0.33 
0 . 36 
0.40 
0 . 45 
0.05 
0.08 
0.13 
0.15 
0.18 
0.22 
0.25 
0.28 
0.32 
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0.35 
0.38 
0.05 
0.08 
0.1 1 
0.145 
0.17 
0 . 20 
0.23 
0 . 26 
0.29 
0.33 
0.35 
0.05 
0.07 
0 . 1 1 
0.14 
0.17 
0 . 195 
0.22 
0 . 25 
0.325 
0.30 
0.33 
0.05 
0.08 
0.11 
0.145 
0 . 17 
0 . 20 
0.23 
0.26 
0.29 
0.33 
0.35 
0.05 
0.03 
0 . 13 
0.15 
0 . 18 
0.22 
0 . 25 
0.28 
0 . 32 
0.3 5 
0.38 
0.05 
0.095 
0 . 14 
0.17 
0.21 
0.25 
0.28 
0.33 
0.36 
0.40 
0.45 
0.05 
0 . 10 
0.15 
0 . 195 
0.245 
0.29 
0.34 
0.38 
0.44 
0.48 
0 . 54 
0.05 
0 . 1 1 
0.17 
0.23 



124 



0 . 29 
0.35 
0.40 
0.46 
0.52 
0.58 
0.64 

//PICT . SYS IN CD * 

RESPONSES WITH HELM STEP COMMAND 
RIBEIRO 

0.0 300.0 - 500.0 300.0 



5.0 



6.0 
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